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Abstract 

We consider an asymmetric control of platoons of identical vehicles with nearest-neighbor interaction. Recent results 
show that if the vehicle uses different asymmetries for position and velocity errors, the platoon has a short transient and 
low overshoots. In this paper we investigate the properties of vehicles with friction. To achieve consensus, an integral 
part is added to the controller, making the vehicle a third-order system. We show that the parameters can be chosen 
so that the platoon behaves as a wave equation with different wave velocities. Simulations suggest that our system has 
a better performance than other nearest-neighbor scenarios. Moreover, an optimization-based procedure is used to find 
the controller properties. 
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1. Introduction 

Control of vehicular platoons has become a field of in¬ 
tensive research. The reason for this is the possibility to 
increase throughput and safety of the highway traffic at 
the same time. Moreover, since the study of such large 
systems gives asymptotic behavior and achievable limits, 
it is also appealing from a theoretical perspective. 

It is well known that if some centralized information is 
present in the system, then the performance is good and 
the system is scalable. These approaches comprise either 
LQR control [T] or local control with added knowledge of 
the desired velocity mu]- Among the systems with per¬ 
manent communication the Cooperative Automatic Cruise 
Control is the most widely studied [4] and implemented. 

In this paper we consider a nearest-neighbor interaction 
without any centralized information. The most important 
performance measure is, apart from the settling time, a so 
called string stability. The system is string stable if the 
disturbance acting at one vehicle does not amplify as it 
propagates along the platoon [ 3 . It is shown in [5] that a 
time-headway spacing policy can achieve string stability. 
However, this policy increases the platoon length with the 
speed of the platoon. 
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A common fixed-distance scenario is a symmetric bidi¬ 
rectional control, where the vehicle measures the distance 
to both its neighbors in the platoon. The properties of 
such control are investigated in HEllH]. One of the ad¬ 
vantages is linear scaling of the T-Lao norm with the number 
of vehicles [8] . The main drawback is a very long transient 
for a higher number of vehicles and sensitivity to noise 
[5]. Both can be qualitatively decreased using a feedback 
control of the leader m- 

Asymmetric control was proposed to shorten the tran¬ 
sient and increase the controllability of large platoons O 
m- a vehicle uses different weights to its two neighbors. 
Nevertheless, the price to pay is an exponential growth of 
Hoo norm of the transfer functions with the graph distance 
between vehicles [I1II31, so this approach is not scalable. 

Both symmetric and asymmetric control share one com¬ 
mon property — the asymmetry in coupling between be¬ 
tween vehicles is the same for position and velocity. That 
is, only one graph Laplacian is used. Except for a simpler 
analysis, there are not many reasons to limit ourselves to 
use only one type of coupling. As was numerically shown in 
|14] and more thoroughly discussed in m, the response 
of a platoon with a symmetric coupling in position and 
asymmetric in velocity can be scalable and have short tran¬ 
sients. Thus, this approach combines the advantages both 
from symmetric and asymmetric control. 

The paper [15] considers a double integrator model of 
the system and uses wave properties for system analysis. 
In this paper we extend the results of [11] to a third-order 
system. The main reason is the presence of friction in ev¬ 
ery real vehicle. We added an integrator to the system to 
achieve coherent solutions. With such a system, the stabil¬ 
ity and performance analysis of the platoon complicates. 
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Although a similar type of a system was recently analyzed 
in m, the results are not applicable for our case since the 
coupling used in our work is non-symmetric. 

Finally, to achieve a good performance, we provide an 
optimization procedure based on the properties of waves 
in the platoon. This approach allows us to optimize the 
controller parameters and weights of the communication 
graph at the same time. 

2. Model of the vehicles 

We assume N + 1 identical vehicles travelling on a line, 
indexed as 0,...,N. The first vehicle with index 0 is a 
leader which is driven independently of the the rest of 
the formation. Unlike standard double integrator mod¬ 
els lanaiis], real systems have a friction, i. e., there is a 
feedback from velocity, which eventually makes the vehicle 
to stop. The vehicle model is 

Xi = —avi + Ui, ( 1 ) 

where Xi is the position of the ith vehicle, Vi = Xi is its 
velocity, a G R is the viscous friction coefficient and Ui is 
the input to the vehicle. 

In order to enable the vehicles in the platoon to track 
the leader moving with constant velocity, we need to sat¬ 
isfy the Internal Model Principle [m ng which in our 
case means the presence of two integrators in the open- 
loop model of each vehicle. Since one integrator is already 
present in the vehicle model (xi = Vi), it suffices to add 
an integral action in the controller of each vehicle. The 
controller is given as q = with defined in ([^ and 
is the state of the integrator in the controller. The input 
to the vehicle is then Ui = Ci. 

Each vehicle uses only the information obtained from 
its nearest neighbors — the vehicle in front of it and be¬ 
hind. The goal of the vehicle is to keep a prescribed spac¬ 
ing to them, i.e., Xi-i — Xi ^ di-iy with dij being the 
desired distance between i and j. The regulation error 
comes from the relative spacing and velocity errors as 

Cz — Qx [(1 di—\^i^ Px{p^i (2) 

-di,»+i)] + 9v [(1 - Pi,)(U-i - Vi) - Pv(vi - Uj+i)] , 

where the position asymmetry is labeled as velocity 
asymmetry as p„ and the gx,gv G R are weights of position 
and velocity errors. The coupling in position is symmetric 
if Pa; = 0.5 and asymmetric otherwise (the same for py). 

To simplify the analysis, we introduce error variables 
Zi = Xi—Xg+do^i- This implies — = Xi-i—Xi — di-i^i 

and ii_i — i* = Xi-i—Xi- We impose that di^^ — dij = dj^k 
and dip = —djp. The single vehicle model combined with 
the controller then has a form 

Zt = -azi + c, Ci = ei. ( 3 ) 

We use a minor state transformation Ci = Ci — aVi to obtain 
a controller-canonical form of the individual-vehicle model 



In a vector form we write the overall system of -I- I 
vehicles (including the leader) as 


z 
z 
z 

(5) 

where z = [zq, ..., zpf]'^. Let us call the system (O a path 
system, since the communication topology is a weighted 
path graph. The Laplacians L,y,Ly G ]^Af-i-ixAr-i-i 
path graph are defined as 
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( 6 ) 
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\ 0 0 0 ... -I 1/ 

/ 0 0 0 0 ... 0 \ 

-(1 - p„) I -Pv 0 ... 0 

0 -(1 - P„) 1 -Py ... 0 


V 0 


0 0 ...- 11 / 


The last vehicle has no follower, so it uses only front 
spacing and velocity errors. This type of boundary condi¬ 
tion is called regular boundary condition [ig . The second 
boundary condition is that the leader is driven indepen¬ 
dently of the platoon (zeros in the first rows of Lj,,Ly). 

We assume that initially the system in is at stand¬ 
still and then the leader starts to move with unit velocity: 

Zi(t) = 0 for i = 1,..., TV, t < 0, (8) 

xo(t) =0, t < 0, xo{t) = t,t>0. 


3. Relation to previous works 


Our paper builds on the results of works [iniiTn]. Both 
papers deal with a signal propagation in systems with 
nearest-neighbor interaction. The vehicle model is a dou¬ 
ble integrator, i.e., Xi = ti with given in ([^, so the 
papers use nearest-neighbor asymmetric interaction. 

The work |19j analyzes a system with a circular topol¬ 
ogy. We call such system a circular system. The interac¬ 
tion between the leader and the vehicle N is added, so the 
Laplacians Lx,Ly G of a circular graph are 


/I -Px 0 

(1 Px) 1 Px 


Px)\ 

0 


.(9) 


V -Px 0 ... -il-px,) 1 / 


Ly is defined similarly by replacing p^, by p„ in ([^. 

For such a system, first the conditions for asymptotic 
stability are found. The most important condition is that 
Px = 0.5 m Prop. 3.5]—there must be a symmetric 
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coupling in the position. For stable circular systems it 
is shown in [m Thm. 4.8] that an external input or a 
disturbance causes two signals to propagate in the system 
in opposite directions and with different velocities. These 
so-called signal velocities are calculated from the phase ve¬ 
locities m Lem. 4.4]. We will use the same ideas in this 
paper (Sec. 4.2 1 to describe a stable system in terms of 
traveling waves. 

The paper m studies transients in path systems. Its 
main result is the description of the transient in the path 
graph using two travelling waves, attenuated at the bound¬ 
aries [m Thm. 3.5]. The connection of the circular system 
from [T5] to the path system (which is the one we are re¬ 
ally interested in) was conjectured in m as follows: the 
asymptotic instability of the circular system should imply 
either flock or asymptotic instability of the path system. 
We now restate flock stability for our system. 


Definition 1 (Flock stability, [l5]). The system is 
called flock stable if it is asymptotically stable and if 
maxtgR \zo{t) — ZN{t)\ grows sub-exponentially in N for the 
conditions given in 


3.1. Assumptions for solving 

The solution in m is based on two main conjectures 
relating the path and circular systems. Although we can¬ 
not prove them, we use them in the present paper as well. 
The final justification of both conjectures is in the match 
of the predicted and simulated values, as shown in Sec. 
The first one states that a local behavior of both systems 
is identical. 

Assumption 1. If the path formation ^ is stable and 
flock stable, then the behavior of a circular system is the 
same as in the path system for vehicles reasonably far from 
the boundaries. 

This assumption allows to use properties derived for 
the circular graph (which is much easier to analyze) in the 
path graph. We note that this treatment implies that the 
boundary condition for the iVth agent in the path graph, 
if reasonable, does not enter the analysis, and therefore 
that boundary condition does not affect our conclusions. 

A similar assumption has been made by others (see 
[DEI [ 20 ]) to simplify the analysis and make the system 
spatially invariant. In fact, in Solid State Physics this idea 
is known as periodic boundary conditions and goes back to 
the beginning of the 20 th century (see ED)- 

In order to investigate flock stability in the path system 
using properties of a circular system, we need an additional 
assumption about relations of the two interconnections. 



(a) Circular system. (b) Path system. 


Figure I: Signal propagation in initial-condition response. 

a) response of an asymptotically unstable circular system, 

b) response of a flock unstable path system. In both cases 

A = 70, a = 3, ga; = 2, = 3, = pa, = 0.33. 


an additional one based on the travelling wave concept 
in distributed control [l0|. Asymptotic instability can be 
caused by the travelling wave which is amplified as it trav¬ 
els in the formation. The amplification will happen far 
from boundaries also in the path system. This results ei¬ 
ther in an asymptotic instability or in a flock instability 
(if the reflections at the boundaries attenuate the signal 
sufficiently) of the path system. 

Both assumptions are illustrated for our third-order 
system in Fig. which shows the initial responses ( 0 : 35 ( 0 ) = 
2 and dij = 0 Vz, j) of the circular and path system. As can 
be seen, the signal gets amplified as it propagates from one 
agent to the other. On the other hand, individual agent’s 
response goes to zero, until the amplified travelling wave 
gets back to the agent after propagating through all other 
agents (see the sharp growth at time 70 in Fig. la I. The 
initial behavior of the path graph and circular systems is 
the same — both amplify the signal. 


4. Analysis of the circular system 

If px ^ Pv in the path system in ([^ , there are two dif¬ 
ferent Laplacians Lx and Ly which are not simultaneously 
diagonalizable. This prevents many convenient approaches 
to guarantee stability such as a synchronization region [ 22 | 
or LMI-based criterion [^. This makes the stability and 
performance analysis of the path system very difficult. 

To overcome this limitation, we invoke the Assump¬ 
tions and to extract some properties of the circular 
system and apply them in the analysis of the path system. 
To obtain such properties, we assume in this section that 
the communication structure is the circular graph with 
Laplacians Lx, Ly. These Laplacians are circulant matri¬ 
ces, which are simultaneously diagonalizable. Note that 
we investigate the circular system only in order to learn 
something about the path system - the circular system is 
not of interest by itself. 


Assumption 2. If the circular formation is asymptoti¬ 
cally unstable, then the path formation is either asymptot¬ 
ically unstable or flock unstable. 

The explanation in m Def. 3.2] is that the path sys¬ 
tem has non-normal eigenspaces which makes the bad ef¬ 
fects more pronounced. To this explanation we can add 


4 . 1 . Stability of the circular system 

When we assume the circular interaction topology, the 
state-space model has a form 
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The Laplacians Lx,Ly are diagonalizable by the dis¬ 
crete Fourier transform. So let Wm be the m-th eigenvec¬ 
tor of L, that is the vector whose j-th component satisfies 

= j = 0,1,..., fV, (11) 

with (j) = 2-Km/{N + \) and i = By [TS], we calculate 

the eigenvalues A^; of gxLx and \y of as 


= gx[^-cos(l) + i{l-2px)siTL(l)\, 

Xv{(t>) = gt,[l — cos())-|-^(l — 2 p„) sini))]. 


( 12 ) 


Let us denote fix = ^ — 2px, /3„ = 1 — 2p„. We can 
expand the eigenvalues Xx and Xy in the Taylor series 

Xxiff}) = gx • ■ ■ ], (13) 

A„(<(>) = gv [^|3v(|) + ■■■]■ (14) 

We now calculate three eigenvalues Vrn,i, * = 1)2,3 
of Mat associated with Wm for each m. The eigenvalue 
equation is 


0 I 0 \ /Um,l 

0 ^ 0 ^ ^ “"*’2 I = V, 

\ gxL/x gv^v dlJ 


m,i I ^m,2 


j 


+ av^ = 0 , 

-I- av^ -I- 2gxV -|- 2gx = 0. 


(18) 

(19) 


The equation (19) implies via Routh-Hurwitz criterion a 


simple necessary conditions for stability. 

Lemma 2. The necessary conditions for the stability of 
(10|) for all N are a > 0, gx > 0 and (?„ > 0 and a > gx/gv 


An important result is that a symmetric interaction in 
position is necessary for the stability of the circular system. 

Lemma 3. The necessary condition for the stability of 
nW for all N is px = 1 / 2 , i.e., Px = 0 . 


Proof. From pT| ) for each value of p we obtain 3 roots 
Two of the roots in (18) are zero. We will find the 
behavior of nearby roots (for |())| small). Since roots of 
polynomials are continuous functions of their coefficients, 
there will be 2 branches of roots with i/(^) tending to zero 
as (j) tends to zero from above, and another 2 branches 
with vffi) tending to zero as (j) tends to zero from below. 

Assume that Px 7 ^ 0. Expand 0 in terms of lowest 
order (in We obtain + iPxPxP = 0. It is easily 

seen that for all values \(j)\ this has two solutions in the 
right half-plane. It follows that if | 0 | is small enough, ( [T7| 
must have two solutions in the right half-plane. It follows 
that we must choose /3a; = 0 , so = 1 / 2 . □ 

The next theorem guarantees the stability of an arbi¬ 
trary large system with a circular topology. 


Theorem 4. All non-trivial eigenvalues of (10) have neg¬ 
ative real part if and only if all of the below hold: 


I. : 

II. : 

III. 


a > 0 and gx > 0 and > 0 and a > gx/gv 
Px = 1/2, 

agv - gx agy - gx\ 


1 - 2 p„ e -- 




Proof. Let us call the statement “Circular system (10) is 


(15) 

from which follows that Um ,2 = Vm,iUm,l, Mm ,3 = Vm,iUm ,2 
and gx^x'^m.l gv^v^m,i'^m,l ~ 

Since Lx and Ly have as an eigenvector Wm, it follows 
from the last equality that itm,i = Wm- For simplicity of 
notation we drop the subscripts of v except from when am¬ 
biguity seems possible. From the last row it is elementary 
to get a set of + 1 eigenvalue equations for Mjq : 

-\- av^ Xyfpfv -\- Xxip) = 0. (16) 

Substituting the expressions for Xx{4>) and Xy{</>), we get 

-\- av^ gy[\ — cosp z(l — 2 p„) sini/Jjz 

-|-pa;[l — cos(/)-|-^(l — 2pa;) sini/] = 0. (17) 

By letting p equal 0 or tt we get real polynomials 


stable” as S. The necessity of conditions I. and II. were 
proved above. We will use them to prove that given I. and 
IT, then III. is false is equivalent to S is false. 

We know from (18) that —a is a solution of 0 and 
that it lies in the left half-plane. By continuity of roots of 
polynomials all the solutions of 0 must lie on a curve 
I'fp) starting at —a. To have unstable roots on the curve 
v{p), the curve must cross the imaginary axis for some p S 
(0, 27r). Then there must be purely imaginary solutions lui 
(w real) to 0. Substitute uv for v into 0 to get 

— auP' [j(l—cos^)—/3« smp]gyUj -I- gx{l—cosp) = 0 . 

( 20 ) 


The real and imaginary parts of (20) are, respectively: 


—— gvPyOJ sin p-\- gx{l — cos p) = 0 , 
— gy{l — cosp)) = 0 . 


( 21 ) 

( 22 ) 


So S is false if both of the last equations hold (i.e., (20) 
has a solution). The equation (22) holds for a; = 0 or 


= 5«(1 ~ cosp), where w = 0 gives only the ‘trivial’ 
eigenvalue (namely p = 0). Plugging the other solution 
= 9v{i — cosp) into ( 21 ) gives: 


/3„ = ± 


ag, 


- gx \/2{l - cosp) 


(23) 


sin/. 

The factor maps the unit circle onto [—oo, — 1]U 

[1, oo]. So for \Py\ > there exists p for which equa- 

tion ( 20 ) is satisfied, the system then has purely imagi¬ 
nary roots and therefore the system can be unstable. If 
\Py\ < ^ then no imaginary solution exists and 

whole curve v{p) lies in the stable half-plane. □ 
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Figure 2: Phase velocities calculated by (24) in vehicles/s 
(solid) and the corresponding damping (dashed) as a func¬ 
tion of (j). There are three waves for each (p, two with 
negative velocity (blue and green) and one with positive 
(red). The phase velocities in green have all high damp¬ 
ing, so they do not affect the signal velocity. The sig¬ 
nal velocities from (26) are shown by red (c+) and blue 
(c_) crosses, the corresponding damping by a red circle. 
N = 500, a = 2,g^ = 6.2, = 10, = 0.4. 


4-2. Signal properties 

Similarly to [T5], we would like to obtain the signal 


velocity in our circular system (10). By [121 Lem. 4.4] 


the phase velocity Cm,i and its damping am,i for a mode 
associated with a given (p = 2Trm/ {N+1) can be calculated 
as 


— Icn(^'y,2,z)/0; ^ 3. (24) 


We are interested in modes with very low damping, since 
they travel in the system with slow decay — they give us 
the signal velocity. Thus, we want to find the eigenvalues 
v{p) with small real parts. They are those corresponding 
to (/i —> 0. To find them, first expand v{p) as 


1 z 

v{p) = imp + -n 2 p‘^ -I- -n^p^ . 


(25) 


We substitute the expansions (131, (141 and (25) into (16). 
Notice that the expansi on ([25 ) works because the terms de¬ 
pending on p cancel in (17). We collect terms of order 0^, 
and p^, etc. The coefficients of these orders must be zero 
and that will determine m. The first non-trivial equation is 
the coefficient of It reads: 0{2) : an\+gyj3yni — 

-9vPv ± \/‘2agx + glPl 


\9x = 


0. We calculate m as m = 


Since for p small by (25) and (24) lin.lv[p)) « mP, the 
coefficient —m determines the signal velocities. 


2a 


Lemma 5. The signal velocities are given as 


c± = 


gy/3y ± \/9vPv + ‘2agx 
2a 


(26) 


where c+ > 0 and c_ < 0 (velocity in vehicles/second). 


By the stability conditions we mentioned, (26) gives 


one positive real and one negative real solutions (red and 
blue crosses in Fig. [^. The wave with the positive velocity 
c+ propagates in the direction with growing vehicle index 
and the wave with c_ the other way. 



time 

Figure 3: Spacing error to the leader zq — Zi with opti¬ 
mized controller (see Sec. [^and|^. N = 250, o = 2,goc = 
6.2, gy = 10 and py = 0.4. 


5. Transients in the path system 


We have obtained enough properties of the circular sys¬ 
tem to derive the transients the original path system (§. 
The transient we analyze is when the platoon is in steady- 
state and the leader starts to move with unit velocity 
We have Theorem [^ guaranteeing the stability of the 
circular system which by Assumption allows for flock 
stability of the path system. The signal velocity in (26) 
should remain the same in the path system — Assumption 
The boundary conditions are the same as in |15j — the 
leader driven independently of the platoon and the agent 
N having no follower. For stable systems (in both senses), 
the orbit of the last agent can be characterized by the fol¬ 
lowing quantities (see Fig. |^: half-period T is the smallest 
t > 0 such that z^ (t) — Zq {t) = 0 and the amplitude Ai of 
the ith oscillation is Ai = max(g[(j_i) 7 ’jzo — zn\. We 
can now restate [13 Thm. 3.5] for the system with friction. 

Theorem 6. If Assumptions [7] and hold and the path 
system is asymptotically stable, the parameter values 
satisfy the conditions in Theorem and as N tends to 
infinity, the system © will behave as a wave equation with 
boundary conditions. The signal velocities are given by 


(26). In particular, if from an equilibrium position at rest. 


the leader starts to move with a unit velocity att — 0, then 
the characteristics of the orbit of Zgit) — z^ft) are: 


A ^ 

Al — I-r, 

|c-hI 


JAfc+i/Afcl = Jc-l/lc+l, 
T = N 


1 1 

|c-h| |c-l 


(27) 

(28) 
(29) 


The proof of [151 Thm. 3.5] uses only boundary con¬ 
ditions and wave velocities, hence it remains valid for our 
case as well. Note that when the leader starts moving, this 
causes first a wave with velocity c+, which then reflects at 
agent N as a wave with velocity c_. Notice that we want 
Jc-l/Jc+l to be less than 1 to avoid exponential growth of 
the amplitudes. Since and a must be positive, we want 
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to keep /3„ > 0, i.e., py < 1/2 and the agent pays more 
attention to the front velocity error. 

6. Simulation verification 


Fig.i numerically validates Theorem by calculat¬ 
ing the relative error between the predicted and measured 
values as a function of N. Let x be a given quantity of 
interest in Thm. § — either Ai, Ai+ijA^ or T. Let Xpred 
be the value predicted by (27), (28) or (29), respectively, 
and Xmeas be the value measured from the numerical sim¬ 
ulations of a finite platoon. The error is calculated as 


^? = log(^^-l) , (30) 

We can see that the relative error of each predicted param¬ 
eter decreases exponentially with the increasing number of 
vehicles in the platoon. This confirms the asymptotic for¬ 
mulas in Theorem and also Assumptions and 



Figure 4: The numerical verification of (27), (28) and (29). 
The relative errors are carried out with (30) for = 6.2, 
Qy = 10, a = 2, px = 0.5 and py = 0.4. The ai and a 2 
are the attenuation coefficients from (28) for A^jA^ and 
A-ijA^., respectively. 


The numerical simulations in Fig. show that a pla¬ 
toon with a controller tuned symmetrically (the left panel) 
has a very long transient. The transient is shortened for 
the case of the asymmetric controller (the middle panel), 
however, the overshoot of such a platoon is extremely 
large, which is a consequence of Lemma — the circular 
system is unstable, hence the path system is flock unstable. 
When we set the asymmetry only in the velocity (right), 
then both the transient and the overshoot are reasonable. 


is based on minimizing the absolute value of the spacing 
error of all vehicles in the formation, denoted as A, when 
the leader starts to move from the stand-still. Therefore, 
the optimization has a form 


^ l‘00 

min £' = min^^ / \ei{t)\dt, (31) 

where the error is given by ei{t) = Zo(t) — Zi(t). Clearly, 
Co = 0. E IS minimized over gx,gv,Pv 

Theorem tells us that the system behaves as a wave 
equation with boundaries. After a unit change of leader’s 
velocity, first the signal spreads from the leader to vehicle 
N with velocity c+ and then it reflects back with velocity 
c_. The graph of the response of the last vehicle in the 
formation must then be almost triangular, as shown in 
Fig.0 The error of the first oscillation for the last vehicle 
(before ejv gets back to zero for the first time) is 


En,i = J \eN{t)\dt « 2 ^^^' 

To get the error of the other agents we assume that the 
maximal value of the error of zth agent is given by (the 
peaks are uniformly spaced from 0 to Ai with distance ^)- 
Then the shape of the error is almost a trapezoid with one 
base of length T and the other with {N—i) = 

T — = T(1—^). The absolute value of the error of 

the ith vehicle in the first oscillation is approximately the 
area of the trapezoid 

(33) 


We have approximated the first oscillations. The errors of 
the others are calculated in a similar same way, i.e. the pe¬ 
riod is again T and the amplitude is obtained using (28). 
The total absolute value of the error of the zth agent is 
Ei = \ei{t)\dt, which is using the trapezoidal approxi¬ 
mation given as the sum of areas of all oscillations 



i 

m 


(34) 


7. Optimization of system parameters 

The previous section gave us signal velocities and am¬ 
plitudes of the transient, which depend on the gains gx and 
gy and the velocity asymmetry py. In this section we give 
an approach how to select these three parameters. We as¬ 
sume that the friction a is given by the vehicle model and 
cannot be affected by the designer. 

We propose the following method for “optimal” (due 
to asymptotic formulas) gain and asymmetry selection. It 


We used (28) to quantify the amplitude of the jth oscilla¬ 
tion and then the sum of geometric series since < I. 


Our criterion (31) captures the sum of Ei of all agents. 


It can be calculated as 


N 


N 


^1 -r 


i=l i=l 

= Ai ^ 


=1 

-TC = JC. 


1 


2n) 1 - ^ 

k+| 


1 _ l£zil 

k+l 


(35) 
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Figure 5: The numerical simulations comparing the responses of three different control strategies for N = 100, when 
the leader changes its velocity from 0 to 1. The figure shows the relative positions to the leader zo{t) — Ziit) of all the 
vehicles for three different combinations of px and p^. For all three cases were px = 6.2, = 10 and a = 2. 


with C = ~ 2 fv) being a constant which cannot 

be changed by optimization. Thus, it suffices to minimize 
J. After plugging from ( [^ and (28), it has a form 

iV2 


J = 


AiT 


|c+| 


|c+| 


1 - V |c+||c-| J |c+| - |c_ 


= JN\{36) 


The number of agents is not part of the optimization and 
does not affect the minimum. Plugging for the signal veloc¬ 


ities from (26) we evaluate the sums and products as |c_ |-|- 
\c+\ = {y/g^pf+2^)/a, \c-\-\c+\ = {gy/3y)/a, |c_||c+| = 
gx/{2a). With these terms the criterion (36) becomes 


J = 


\/5p| 


9vPi 


2a^2a _ M 


2a 


9lPv9'~ 


-2a. 


(37) 



Figure 6: Logarithm of E{N) for different architectures 
— symmetric with Px = Pv = 0.5 (line 1), asymmetric 
with Px = Pv = 0.4 (line 2) and asymmetric with px = 
0.5,= 0.4 (line 3). Line 4 shows estimate using (40). 
Other parameters were px = 6.2, py = 10, a = 2. 


Since a is a given constant and the square root is mono¬ 
tonic function, we get the final optimization problem 


min 


1 2a 
9x g‘iPl9x ’ 

s.t. conditions in Theorem ID 


(38) 


The final criterion is a function only of the parameters 
9x,9v and /3„, which the platoon designer can affect. 


The scaling of E of different architectures with N cal¬ 
culated from simulations is in Fig. It is clear that 
the error is the smallest for asymmetry only in velocity 
and also that the error of asymmetric control with iden¬ 
tical asymmetries scales exponentially in N, which con¬ 
firms flock instability {px yf 0.5). Also the predicted value 
matches the calculated one. 


7.7. Scaling of the absolute error 


The total error (31) can be written using (35) and (36) 

N 


E^JC = JN"^ 


i-1 




N \ 2N J 


(39) 


Lemma 7. The error E in (31) scales cubically with N 
as 

EiN) « + 1){4N - 1). 


(40) 


Proof. The proof is a simple manipulation of (39). 

N 


N 




N 


N 


2^E*-E* 

i=l i=l 

N{N + 1){2N + 1) 
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= 7|2iV 
J 


2N^ 

=1 

N{N+1) 


(41) 


=—N{N+ 1){4:N-1). □ 


8. Optimization results 


Although both the criterion (38) and the stability con¬ 


straints are nonconvex, the optimization using function 
fmincon in Matlab for nonconvex optimization terminated 
quickly and successfully. The code used for simulations in 
the whole paper can be obtained at |M| . 


As follows from (38), the optimization procedure tried 
to increase the gains px and py to decrease the criterion. 
Therefore we specified upper bounds on PxtPv to limit the 
controller effort. The optimization was conducted for a 
given friction a = 2 and we got the values px = 6.2, py = 10 
and Py = 0.4. The upper bounds for both gains px, Py were 
set to 10. To stay away from the flock stability boundary, 
we changed the flock stability criterion in Thm. to 

p 1 r, ^ ( a9v-9x , ^ apy-Px 

j’ ' ’ 

with e = 0.1. The response is shown in Fig. [^and[^. 
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(a) Sens, to change in a. (b) Sens, to change in pt,. 

Figure 7: Sensitivities to different parameters. The plots 
show value of i? of (31) calculated from simulations. 


8.1. Robustness evaluation 

The optimization results should be verified to give ro¬ 
bust results. The simplest way to achieve robustness is to 
add some nonzero term e to each of the stability criteria in 
Thm. 1^ similarly to (421. Then the system is not allowed 
to operate on the flock stability boundary. 

The most important parameter of the system is the 
friction a. This might change during the operation of the 
system and also might not be exactly known apriori. Us¬ 
ing the values gx = 6.2, = 10, p„ = 0.4,7V = 1200, 

we simulated the response of the system for friction range 
a G [1.4, 2.8] and calculated the norm of the error using 
(Hj. Fig. |7a| shows how the norm changes with friction. 
We see that the change is approximately linear in friction 
and the system has a good performance for a wide range 
of a. The sharp growth for low friction caused by a flock 
instability confirms the Assumption]^ The stability crite¬ 
rion of a circular system in Thm. |^is violated for a < 1.514 


since = 1 — 2p„ > 




making the system flock 


unstable. As TV —oo, the sharp growth appears exactly 
at the critical point a = 1.514. 

The Figure [7b| shows how the error norm changes as a 
function of p„ € [0.34,0.44] with the optimal value p* = 
0.37 (other parameters are gx = 8.3, = 10, o = 2,iV = 

1200). These are the values with e = 0 in ( [4^ . It is clear 
that the value p* is almost the minimum of the function 
E{py). The better performance for lower (non-optimal) 
Pv < pt is due to asymptotic formulas used. When number 
of vehicles increases, the minimum will get closer to p*. 
Due to the sharp growth for p„ < 0.36 we recommend 
using p„ = 0.40 to achieve robustness, as obtained using 


(42) with e = 0.1. Then we get also gx = 6.2, as above. 


9. Conclusion 


We investigated the transients in a vehicular platoon 
with identical vehicles having friction. The asymmetries 
for a coupling in position and velocity are different. For 
some steps in the analysis, we considered circular system 
in order to infer properties of the path system. 

To achieve flock stability, the coupling in the position 
must be symmetric. Then the path system behaves as a 
wave equation with boundaries, where the travelling waves 


have different velocities. Using the velocities, we developed 
an optimization procedure to determine the parameters of 
the controller. We conclude that the behavior of a platoon 
with different asymmetries is superior to both symmetric 
and identically asymmetric platoon. 

We believe that the approach and some of the results 
shown here easily generalize to more complicated systems 
having two integrators in the open loop. 
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